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data collected at the T{AS) resonance with the BABAR detector at the PEP-II e'^e" storage rings. 
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mixing. We obtain a 90% confidence level upper limit of B{B -> AtXrVi)/B{B -> At X) < 3.5%. 

PACS numbers: 13.20He, 12.38.Qk, 14.40Nd 



Decays of B mesons to charmed baryons are not as well 
understood as those to charmed mesons. In particular, 
there is limited knowledge, both theoretical and experi- 
mental, about semileptonic B decays to the charmed 
baryon [l|. li B decays to charmed baryons are domi- 
nated by external W emission (Fig. [IJi) , as is the case 
for B decays to charmed mesons [2|, y], and final-state 
hadronic interactions are small, the semileptonic fraction 
of these decays should be roughly the same: 

B{B A+Xi-Ve) B{B DX"t-Vi) 
B{B^ A+X') B{B DX'") ^ ' 

where i = e or fx, and D is understood to be or 
£)(*)+ The semileptonic fraction of B decays to charmed 
mesons is currently measured to be 11.1 ± 0.8% A 
significantly smaller semileptonic ratio for B decays to 
charmed baryons would be evidence for a sizable internal 
W emission amplitude in baryonic B decay (Fig. [T]d) or 
a large contribution of final state interactions. 

About 90% of the measured inclusive semileptonic 
B — )• Xcl^T^i branching fraction into charmed final states 




u M u d 



FIG. 1: Feynman diagrams for B decays into a charmed 
baryon through external W emission (a) and internal W emis- 
sion (b). 



can be accounted for by summing the branching fractions 
from exclusive B — ^ D'^*\-K)l^'Ut decays 5]. Semilep- 
tonic B decays to charmed baryons could account for 
some of the remaining difference. 

A previous search for semileptonic B decays into 
charmed baryons by the CLEO collaboration Q re- 
sulted in an upper limit on the ratio B{B — > 
A+Xe'Ve)/B{B/B A+X) < 5% at the 90% con- 
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fidence level. By using B{B/B A+X) = 0.045 ± 
0.004 ± 0.012 and 3(8 A^X)/B{B A+X) = 
0.19 ± 0.13 ± 0.04 [3], and assuming lepton universality, 
this result implies a semileptonic fraction limit B{B ~^ 
A+Xt-Vt)/B{B ->■ A+X) < 6% at 90% confidence level. 

There are two caveats to the CLEO measurement. 
First, the electron candidate is required to have a mo- 
mentum greater than 0.6 GeV/c, which reduces back- 
ground due to fake and secondary electrons, but may 
also reduce signal efficiency. Second, because the CLEO 
measurement was unable to constrain the flavor of the B 
meson, the quoted fraction suffers from large systematic 
uncertainties due to the uncorrelated B A^ X back- 
ground. We address these two points by reconstructing a 
B meson in a hadronic mode and look for the signal in its 
recoil. The resulting sample has less background, which 
allows us to lower the lepton momentum cutoff, and the 
flavor of the hadronic B meson determines the flavor of 
the signal B, up to mixing effects. By normalizing to 
the correlated B — >■ A+X decay mode, many systematic 
uncertainties cancel. 

In this paper, we present a search for semileptonic B 
decays to A+ using data collected with the BaBAR de- 
tector at the PEP-II asymmetric-energy e+e~ storage 
rings at SLAG. The data consist of a total of 420 fb~^ 
recorded at the T{4S) resonance between 1999 and 2008, 
corresponding to approximately 460 million BB pairs. 
The BaBAR detector is described in detail elsewhere 0- 
Charged particle trajectories are measured by a five-layer 
double-sided silicon vertex tracker (SVT) and a 40-layer 
drift chamber (DCH), both operating in a 1.5 T mag- 
netic field. Charged particle identification is provided by 
the specific ionization energy loss (dE/dx) in the track- 
ing devices and by an internally reflecting ring-imaging 
Cherenkov detector (DIRC). Photons are detected by 
a CsI(Tl) electromagnetic calorimeter (EMC). Muons 
are identified by the instrumented magnetic-flux return 
(IFR). A detailed Geant4-based Monte Carlo (MC) sim- 
ulation Q of BB and continuum events (light quarks 
and r pairs) is used to study the detector response, its 
acceptance, and to test the analysis techniques. 

We search for semileptonic B — A+X£~17i decays with 
£ = e or fi in events pre-selected to contain a candidate B 
reconstructed in a fully hadronic decay mode (i?tag), as 
described later in the text. We select signal candidates 
in these events by looking for candidate leptons and fully 
reconstructed A+ decays. We then refine our selection 
of Stag, and make a final signal extraction based on the 
selected Btag and A+ kinematic properties. We also se- 
lect candidate B A+X events, starting with the same 
sample and using similar techniques and selections, but 
without requiring an identified lepton candidate. 

Selection criteria are optimized using MC simulation 
of signal and background processes. Because little is 
known about B A+Xt~vi decays, we use a signal 



model which can be tuned to cover a large range of pos- 
sible kinematics of the final-state particles. In this model, 
the B decays semileptonically into an intermediate mas- 
sive particle Y , B ^ Y£~17£, with a kinematic distribu- 
tion according to phase space @. The Y subsequently 
decays into a A+, an anti-nucleon (anti-proton or anti- 
neutron), and ni (122) charged (neutral) pious, again as- 
suming phase space distributions. The free parameters 
in the model (the mass my and width Ty of the pseudo- 
particle Y, and ni and 712) are tuned to reproduce the 
lepton and charmed hadron momentum spectra predicted 
by the B L»(*)7r£l7£ model of Goity and Roberts 
after accounting for the phase space limits implied by 
the large baryon masses. We choose my = 4.5 GeV/c^, 
Ty = 0.2 GeV/c^, and ni + 713 < 6. 

We reconstruct Stag decays of the type B — >• DY' , 
where Y' represents a collection of hadrons with a total 
charge of ±1, composed of n[TT+ + n'2K+ + n'-^K'^ + 71,471°, 
where n'^+ n'2 < 5, rig < 2, and 71,4 < 2. Kg candidates 
are reconstructed in the tt+tt" decay mode, tt^ candidates 

Q 5):Q 

in the 77 mode. Using D (D^) and D {D*~) as seeds 
for B+ {B°) decws, we reconstruct about 1000 complete 
B decay chains [lll|. 

The kinematic consistency of a Stag candidate with 
a B meson decay is evaluated using two variables: the 
beam-energy substituted mass uies = \/s/4 — jp^jp, 
and the energy difference AE = Eg — -y/s/2. Here ^/s 
is the total center of mass (CM) energy, and p*g and 
E^ denote the momentum and energy of the Stag candi- 
date in the CM frame. For correctly identified Stag de- 
cays, the rriES distribution peaks at the B meson mass, 
with a resolution of about 2.5 MeV/c^ averaged over the 
decay modes, while AE is consistent with zero, with a 
resolution of about 18 MeV. We select Stag candidates 
in the signal region defined as 5.27 GeV/c^ < tues < 
5.29 GeV/c^, with a AE within 4(t of zero. This selec- 
tion has an estimated efficiency of 0.2% to 0.3% per B 
meson. 

We identify electron and muon candidates by combin- 
ing the information on the measured momentum and en- 
ergy loss in the SVT and DCH, the angle of Cherenkov 
radiation in the DIRC, and the energy deposition and 
shower shape in the EMC. For sufficiently hard muons, 
the information from the IFR is also used. We correct 
for bremsstrahlung of electrons by combining the four- 
momentum of the electron with those of detected pho- 
tons which are emitted close to the electron direction. 
We require lepton candidates to have a momentum in 
the CM frame > 0.35 GeV/c and a point of closest 
approach to the collision axis of less than 0.1 cm. The 
p*i selection value is motivated by the large mass of the 
A+ and the assumption of another baryon in the decay 
due to baryon number conservation, which greatly re- 
stricts the kinetic energy available to the leptons. We 
identify photon conversions and i:^ Dalitz decays using a 
dedicated algorithm based on the vertex and kinematic 
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properties of two opposite charge tracks, and eliminate 
electron candidates coining from these. 

Candidate baryons are reconstructed in the 

pK~'K'^, pKg, pKg'iT~^'!T~, ylTT"*", and yl7r+7r''"7r~ modes. 
A candidates are reconstructed in the pTr~ decay mode. 
Only Ac candidates with opposite charge of the lepton 
candidate are considered. Charged daughters of the 71+ 
candidate are fit to a vertex tree with Kg and A 
masses constrained to their known values Q , and the 71+ 
origin constrained to the known average luminous posi- 
tion of the beams within its measured size and uncertain- 
ties. In events with multiple yl+ and/or £ candidates, the 
candidates are fit to a common vertex, and the A'^£~ pair 
with the highest vertex fit probability is selected. 

We refine the selection of Stag candidates by first 
removing those whose daughter particles are based on 
tracks already used to reconstruct the signal-side 71+ or 
lepton and those charged iJtag candidates whose flavor 
is opposite that of the signal B candidate. We account 
for mixing effects by weighting B'^ and B^ tags accord- 
ing to the Ac charge, as described in Ref. [13|. In events 
with multiple -Stag candidates, we select the one recon- 
structed in the highest purity mode, where the purity is 
estimated for each i?tag decay chain using MC simula- 
tion as the ratio of signal over background events. When 
multiple candidates in the same event have the same i^tag 
mode, we select the one with the smallest \AE\ value. 

We reconstruct the CM missing momentum Pmiss by 
noting that Pmiss + Pvis = in the CM frame, where the 
visible momentum pVis is computed by summing the mo- 
mentum vectors of the Stag, the Ac and £ candidates, plus 
any additional well measured charged track or neutral 
cluster boosted to the CM frame. We require |pmiss| > 0.2 
GeV/c to remove background from hadronic B — ^ tI+X 
decays in which all the particles in the X system have 
been reconstructed and one hadron is misidentified as a 
lepton. We compute the total observed charge of the se- 
lected events by adding the charges of all particles used 
in the Pmiss calculation, and require this to be zero. This 
reduces the background in the i?tag reconstruction due 
to missing particles. 

Backgrounds are divided according to whether they 
contain a correctly reconstructed 71+ candidate. Those 
which contain such a candidate are called "peaking back- 
ground," while those that do not are called "combinato- 
rial background." The predictions from MC simulation 
of generic BB and continuum events show that the peak- 
ing background arises mainly from hadronic B — > 7l+X 
decays, where the 7l+ is correctly reconstructed, and 
the lepton candidate is an electron from gamma con- 
Dalitz decays, or a hadron misidenti- 



versions or tt 
fled as a muon: 



we estimate 3.6 ± 0.7= 



± 0.7, 



syst. 



and 



15.3± 1.5stat. ± l-4syst. peaking background events for the 
electron and muon samples, respectively. The relatively 
large peaking background rate for the muon channel is 
due primarily to the low lepton momentum cut. 



We determine the B semileptonic signal yield with a 
simultaneous unbinned maximum likelihood fit to the dis- 
tribution of the 7l+ invariant mass on both the electron 
and muon samples. The 7l+ invariant-mass distribution 
is described by the sum of three probability density func- 
tions (PDFs) representing signal, peaking background, 
and combinatorial background. The functional forms of 
the PDFs are chosen based on simulation studies. The 
signal and peaking background contributions are mod- 
eled as Gaussian functions whose mean and width are 
fixed to the values obtained from a fit to the the 7l+ 
candidate mass spectrum in the B — !■ 7l+X data sam- 
ple described below. The number of peaking background 
events is fixed to the prediction from MC simulations. 
The combinatorial BB and continuum backgrounds are 
modeled as a first-order polynomial, whose parameters 
are constrained by a fit to the 7l+ invariant mass side- 
bands, defined as the mass ranges from 2.23 — 2.26 and 
2.31 - 2.34 GeV/c^. The fit to the 7l+ invariant mass 
is shown in Fig. [21 projected separately for the electron 
and muon samples. The corresponding yields are shown 
in Table m 
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FIG. 2: Fit to the 7l+ candidate mass distribution for B — ^ 
A'^Xe~Ve (a) and B — >■ A^ X fj,~Vfj, (b). The data are shown 
as points with error bars, the overall fit as a solid line, and 
the peaking background contribution as a cross-hatched area. 
The combinatorial BB and continuum background is shown 
as the area below the dotted line. 

In order to reduce systematic uncertainties due to -Btag 
and 7l+ reconstruction, the B — 7l+X£~T7£ branching 
fraction is measured relative to the inclusive B{B — 
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A'^X) branching fraction. To determine the inclusive 
yield, we start with the same i^tag sample used for the 
semileptonic selection. We reconstruct yl+ candidates as 
in the semileptonic case, choosing the candidate with the 
highest vertex probability in case of multiple candidates. 
We exclude -Btag candidates with daughter particles in 
common with the A'^ candidate and resolve multiple -Btag 
candidates as in the semileptonic case. 

We determine the B — >■ A'^X signal yield with an un- 
binned maximum likelihood fit to the A'^ invariant mass. 
The fit function consists of the sum of two PDFs repre- 
senting signal and combinatorial background, described 
by a single Gaussian and a first order polynomial, respec- 
tively. All parameters of the signal Gaussian are left free 
in the fit. We obtain a A+ mass value of 2.2853 ± 0.0003 
GeV/c^, consistent with the current world average 
and a resolution of 4.0 ±0.3 MeV/c^, consistent with ex- 
pectations from MC simulations. The A'^ invariant mass 
distribution on the inclusive sample and the results of the 
fit are shown in Fig. [31 




2.24 2.26 2.28 2.3 2.32 2.34 

Ac candidate mass (GeV/c^) 



FIG. 3: Fit to the candidate mass distribution for B — > 
A+X. The data are shown as points with error bars, the over- 
all fit as a solid line, and the combinatorial BB and continuum 
background as a dashed line. 

We determine the relative branching fraction B{B — > 
A'^X£~Vi)/B{B — )• A'^X) as the ratio of the measured 
signal yields, after correcting for the ratio of the recon- 
struction efficiencies: 

B{B^A+Xn) ^ /NA feA 

B{B^A+X) \Nj\eJ- ^' 

Here, (Ni) is the number of B ^ A+Xi-Vi (E 
A'^X) events reported in Table |T] together with the cor- 
responding reconstruction efficiencies Cs (ei); the latter 
include the i?tag efficiencies, which are estimated with 
MC simulation. 

Many systematic uncertainties approximately cancel in 
this ratio, such as those due to the yl+ and i?tag recon- 
struction efficiencies and the A'^ decay branching frac- 
tions. We categorize the remaining systematic uncertain- 
ties into those which directly affect the signal yield, and 



those which affect only the efficiency. The systematic un- 
certainties that have been considered are described below 
and summarized in Table ITll 

Systematic uncertainties in the signal yield are dom- 
inated by the peaking background yield. We estimate 
this by propagating the uncertainty in the B — ^ A'^X 
branching fraction, and the Poisson error from the MC 
simulation. We add in quadrature the effect of varying 
the probability for a pion to be misidentified as an elec- 
tron or as a muon by 15%, where the range is estimated 
using data control samples [lH . Systematic uncertainties 
due to background electrons from photon conversions and 
tt" Dalitz decays are negligible. 

To account for a possible bias due to the fit technique, 
we prepare ensembles of MC experiments, in which events 
are generated according to the PDF shapes determined 
from data. We vary the signal to background ratio and 
fit for the signal as in the full analysis. The average 
difference between the fitted value of the yield and the 
true value is taken as a systematic uncertainty, labeled 
"Fit bias" in Table HI 

TABLE I: Signal yields and reconstruction efficiencies for the 
B AtXrVe, B AfX, and B/B -^• yl+X decays with 
the corresponding statistical uncertainties. 

Decay mode iVdata e (xlO~^) 

B AtXe'Ve 15.0 ±6.8 1.98 ± 0.17 
B ^ A+Xii-Uf, -6.2 ±6.3 1.04 ±0.12 
B^A+X 934 ± 55 3.09 ±0.11 

B/B ^ A+X 1386 ±66 3.21 ±0.12 



Systematic uncertainties on the reconstruction effi- 
ciency ratio are dominated by the uncertainty in the sig- 
nal model. This is estimated by comparing our nominal 
signal model with a pure phase space model, where the 
B — )• A+Xi~V^ decay occurs in one step, taking the full 
difference in the signal efficiency estimate compared to 
our nominal signal model as the systematic uncertainty. 
The larger systematic uncertainty for the muon channel 
is due to the low muon identification efficiency for the 
soft leptons. The uncertainty in the reconstruction effi- 
ciency due to the limited statistics of the MC simulation 
is added as a systematic uncertainty by weighting the 
events to the data size. The peaking background in the 
inclusive mode due to cc is estimated using the prediction 
from our MC simulation and is found to be compatible 
with the statistical uncertainty of the sample, which we 
take as a systematic uncertainty. We estimate the sys- 
tematic uncertainty on the signal efficiency due to par- 
ticle identification by varying the electron (muon) iden- 
tification efficiency by 2% (3%), based on studies using 
data control samples [nj. Since the order for selecting 
the best candidate is different between the semileptonic 
and inclusive samples, the uncertainties on the ratio of 
the Stag and A+ efficiencies do not exactly cancel. We 
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evaluate the corresponding systematic uncertainty by re- 
versing the order of the lepton and -Stag selection and 
comparing with our standard selection order using the 
same MC simulation of our signal model used to estimate 
the reconstruction efhciency. Since we find the reversed 
selection order efficiency to be compatible with the stan- 
dard selection order efficiency within the precision of our 
MC simulation, we estimate the systematic uncertainty 
as the statistical uncertainty of that comparison. 



TABLE II: Sources of systematic uncertainties. 



Yield systematics (events) 




£ = 


Peaking background: sample statistics 


1.0 


1.4 


Peaking background: B{B J^t^) 


1.6 


4.7 


Lepton mis-id rate 


0.7 


2.0 


Fit bias 


0.3 


1.2 


Total 


2.0 


5.4 


Efficiency ratio systematics (%) 


i^e 


£ = fi 


Signal model 


11.3 


35.9 


Reco. efficiency statistics 


8.4 


11.4 


Peaking background: B — >■ A'^X 


1.9 


1.9 


Lepton id efficiency 


1.1 


2.7 


Selection order 


5.0 


6.8 


Total 


15.1 


38.4 



The central values of the branching fraction ratios are 
summarized in Table IIIIl We find a signal significance 
S = 2.1, including the systematic uncertainties on the 
signal yields, from the difference in the log likelihood val- 
ues between the nominal fit and a fit in which we fix the 
signal yield to zero. By scanning the likelihood values 
including the full systematic uncertainties, we estimate 
an upper limit at the 90% confidence: 



B{B A+Xl-Vj) 
B{B^A+X) 



< 3.5%. 



(3) 



For a comparison with the CLEO result 6], in 
which the flavor of the semileptonic B was not de- 
termined, we repeat the analysis without requiring 
the charge-flavor correlation between the i?tag and the 
yl^ in the inclusive mode. The corresponding yield 
for the inclusive mode is shown in the last row of 
Table ID We obtain the branching fraction ratio 
BiB ^ A+Xi-ue)/BiB/B ^ A+X) = (1.2 ± 0.7stat. ± 
0-4syst.)% with its corresponding 90% confidence level up- 
per limit B{B A+Xe-Vi)/B{B/B A+X) < 2.5%, 



TABLE III: Central values of the branching fraction ratio 
B{B At Xrvi)/B(B A+X). The last line averages over 
e and /i. 



Mode 


BE ratio 


l = e 


2.5=bl.lstat. i0.6syst. 


i = fi 


■-2.0±2.0stat.±1.9sy.t. 


I = e, fi 


1.7±1.0stat.±0.6.y.t. 



which improves the CLEO limit. We find that remov- 
ing the charge-flavor correlation between the lepton and 
the Stag in the semileptonic mode also yields consistent 
results after re-estimating backgrounds. 

In conclusion, we have presented a search for semilep- 
tonic B decays into the charmed baryon A+. We ob- 
tain an improved upper limit with respect to previ- 
ous measurements [_6I| on the relative branching fraction 
B{B ->■ A+X£-ue)/B{B A+X), which is found to be 
much smaller than the corresponding relative branching 
fraction for B decays into charmed mesons. Our result 
shows that the rate of baryonic semileptonic B decay is 
too small to contribute substantially to the branching 
fractions of inclusive semileptonic B decays. 
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